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1 Introduction 

The quadratic form of the increments of a diffusion type process under finite time horizon, the "realized volatil- 
ity" in financial context for example, is in general asymptotically mixed normal. When the limit is normal, 
the asymptotic expansion of the quasi-likelihood type estimator was derived in [6] as an application of the 
martingale expansion. The expansion for the quadratic form with asymptotically mixed normal limit is then 
indispensable to develop the higher-order approximation and inference for the volatility. However, the classical 
approaches in limit theorems, where the limit is a process with independent increments, do not work. 

The previous paper [8]^ presented asymptotic expansion of the martingale with mixed normal limit. The 
expansion formula is expressed by the adjoint of a random symbol with coefficients described by the Malliavin 
calculus, differently from the standard invariance principle. As an application, an asymptotic expansion for a 
quadratic form of a diffusion process was derived in [8] . The aim of this article is to give a short review of this 
result and some details of its derivation. 



2 Asymptotic expansion of a quadratic form of a diffusion process 

We consider a diffusion process satisfying the Ito stochastic integral equation 

Xt = Xo+ f b{Xs)ds+ f a{Xs)dws. (1) 

Jo JQ 

Here b and a are assumed to be smooth functions with bounded derivatives of positive order. We only treat 
one-dimensional case for notational simplicity, however, multivariate analogue is straightforward. Extension to 
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Sciences, University of Tokyo; and by a Cooperative Research Program of the Institute of Statistical Mathematics. 
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general Ito processes is also possible but the deseriptions would become involved. We will consider a quadratic 
form 



(2) 



of the increments of X, where tj = j/n and AjX = Xt- — Xt-_i . The function c is in Cj?°(M).^ 

The quadratic form (2) of the increments of X appears in applications in statistics and finance. In the 
high-frequency sampling of n tending to oo, Un converges in probability to 



Uoo = \ c{Xt)a(Xtfdt. 



The normalized error is 



(3) 



It is well known that Z„ has a mixed normal limit distribution in general. However, the limit theorem is not 
always sufficient for approximation nor for theoretical statistics. Our interest is in more precise approximation 
to the distribution of Z„. 

We write ft for f{Xt), given function /. For differentiable /, the Ito decomposition of ft = f{Xt) is denoted 

by 



ft — ./o 



ffUs. 



Obviously, 



ff^ =a{Xt)d,f{Xt) and ff^ = Lf{Xt) with L = 69, 



For a cZi-dimensional reference variable, we will consider 

1 " 

-F„ = - or Fr,=F^ 



(4) 



where /3 e Cf'(R;R''i). The results will be the same in these cases up to the first order asymptotic expansion 
we will discuss in this paper. It is standard in theoretical statistics to treat the joint distribution of Z„ and 
because F„ can be the quadratic variation of the score martingale and the LAMN property is then established 
on the joint convergence. The studentization also motivates the joint convergence. 
Let a{x) = c{x)a{xY . Define 1 + di-dimensional vector fields Vb and Vi by 



b{xi) - \(j{xi)dxi(T{xi) 



and Vi(xi,X2) 



a[xi) 




for xi G R and X2 G R''^. Denote by Lie[Vb; Vi\{xi^X2) the Lie algebra generated by 

Vi, [V,,Vj] (*,j=0,l), [V,,[V,,Vk]] (z,j,fc = 0,l),.... 

at {xi^X2), where [•, •] is the Lie bracket. 

Assume that the support suppP"^" of the law of X^ is compact. We will assume the following non degeneracy 
conditions. 



[HI] 



inf,. 



\a{x)\ > 0. 



is simply denoted by 



) is the set of M'^-valued smooth functions defined on R'' with all derivatives of at most polynomial growth. 
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[H2] Lic[yo; Vi](Xo,0) = Ri+^'i a.s. 



The asymptotic expansion formula will be described with certain random symbols. The full random symbol, 
denoted by (t(z, iu, iv) for (z, m, u) e M x R x M^^i , consists of the adaptive random symbol a_ and the anticipative 
random symbol a. These random symbols arc defined as follows. 

Let 

u a , 1 [0] 2 [1] [1] 

and 

kt = 2ctht(Jt + ct(Jt(y^t^ - -cf'a-'l. 
In the present situation, the adaptive random symbol is given by 

g_{z,iu,iv) = — / a{Xs)^ds( f a{Xs)^ds\ {iu)^+iu f ktdwt + iu f htdt. (5) 
3 Jo ^Jo ' Jo Jo 

The processes DgXt and DrDgXt satisfy a system of partially linear equations: 

DsXt = a{Xs)+ f b'iXt)DsXt,dh+ f a\Xt,)D,Xt,dwt, 

J s J s 

for s <t^ and 

DrD.Xt = a'{X,)DrX,+ f b"{Xt,)DrXt,D,Xt,dh+ f b'{Xt,)DrD,Xt,dh 

J s J s 

+ ^* a"{Xt, )DrXt, DsXt, dwt, + ^* a\Xt,)DrD,Xt, dwt, 

for r<s < t. The L^-estimates of the solution are at hand. Then the anticipative random symbol is given by 
the formula 



cr(m, iv) 



iu a{Xs)(7s,siiu,iv) ds 



(6) 



with 



(T,,,(^u,w) = (^-u^ I a'{Xt)DMt + i I l3'{Xt)[v]DMt^ 



-u^ f {a'\Xt){D,Xt)^ +a'{Xt)DsD,Xt}dt 

J S 

+1 [ {l3"{Xt)[v]{DsXtf + p'{Xt)[v]DsD,Xt}dt, 



(7) 



where a{x) = a{x)'^. Let Coo — 2 a{Xt)dt. 

The density of the multi-dimensional normal distribution with mean vector m and variance matrix C is 
denoted by m, C). With the full random symbol 



(7{z,iu,iv) = a(z,iu,iv) + a{iu,iv)^ 
the density function pn{z,x) £ C°° {M}^'^'^) is defined by 

1 



Pn{z,x) 



E 



(z;0,Coo)4(Foo) 



-.E 



a(z,a„a,)*U(z;o,Coo)4(Foo) 



(8) 



(9) 
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Here 5x{Foo) is Watanabc's delta function ([5]). The adjoint operation (T{z,dz,dx)* is defined by 



j;0,Coo)4(Foo)| = ^(-a:,n(-9,)"^|c,(z)0(z;O,Coo)<5,(i^oo)| 



for the random symbol cr(z, dz, dx) having a representation 

(j{z,dz,dx) = ^Cj(z)(m)™^(w)"^- (finite sum) 

3 

where Cj are random functions of z, rrij G Z'J^ (in the present case d ~ 1) and rij G Z'^^ . If Cj, Coo and i^oo are 
smooth in Malliavin's sense and F^o satisfies a suitable nondedeneracy condition, Coo being nondegenerate as 
well, then this adjoint operation is well defined. These conditions are satisfied in the present situation, therefore 
Pn{z,x) is well defined. See [8] for details of random symbols and the adjoint operation. 
The following theorem gives an error bound for the approximate density Pn{z,x). 

Theorem 1. Suppose that [HI] and \M2\ are satisfied. Then for any positive numbers M and 7, 



E[fiZ„,F„)] - / f{z,x)pn{z,x)dzdx 



sup 

E£(A/.- 

asn—>oo, where£{M,j) is the set of measurable functions f : B}^"^^ R satisfying \f{z,x)\ < A/(l + |z| + |a;|)''' 
for all {z,x) e K X R''i. 



Remark 1. Condition [HI] is usually ensured by the uniform ellipticity of the diffusion process Xt and a 
reasonable choice of the estimator. In this sense, it is a natural assumption in statistical context. 

Remark 2. The hybrid I method (a rough Monte-Carlo method in the first order asymptotic expansion term) 
is useful in the application of the expansion formula to numerical approximation. Applications to volatility 
derivatives are in our scope. 

Remark 3. In the present article, we have a conditioning variable as F„. On the other hand, it is also possible 
to consider versions of our results without Fn if we were interested in a single (not joint) expansion. It will 
reduce the differentiability conditions of variables. However, considering the joint distribution is natural in 
non-ergodic statistics. Studcntization is important in any case. 

Remark 4. It is also possible to obtain asymptotic expansion of the conditional distribution. 
Section 4 will give some details of derivation of Theorem 1. 

3 Review of the asymptotic expansion of a double stochastic integral 
having a mixed normal limit 

In this section, we will give a short review of the martingale expansion. We refer the reader to [8] for details. 

On a stochastic basis (il, J^, F = {Tt)te[o,i]i P) with T = Ti, we shall consider a sequence of d-dimensional 
functionals with decomposition 

Z„ = A/„ + VF„+r„7V„. (10) 

Here, for every n G N, M" = (Af")(g[Q denotes a d-dimensional continuous martingale with respect to F and 
Mn = M". In this decomposition, we assume Wn, Nn € J^(fi; M."^)^ and (r„)„N is a sequence of positive numbers 
tending to zero as n — > 00. 

In this section, the reference variables are F„ S J^(il;M'*i) {n G N) are general and we do not assume a 
specific structure like (4). It is possible to give asymptotic expansion of £{Z„} under certain conditions; see 

^The set of d-dimensional measurable mappings. 
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[8]. The same paper apphed the expansion to the case where M„ is given as a sum of double Ito integrals, as 
reviewed in what follows. 

Let (W, P) be an r-dimensional Wiener space over time interval [0,1] and let H be the Cameron-Martin 
subspace of W. Suppose that the probability space (fi, J", P) is such that = fi' x W, J" = J"' x B(W) and 
P = i-" X P for some probability space (il', J-"', P'). We will use the partial Malliavin calculus on Vt based on 
the shifts in the direction of H . For a Hilbcrt space E, the Sobolev space of E'-valued functionals on $7 with 
indices s G R for differentiability and p E (1, oo) for integrability is denoted by p(£'). 

Let d = d+di andi = d + 6. Let e Di+i^^{H (g) R'^) and isT" S Bi+i,oo{H (g)W (g) R''). Since H can 
be identified with i^([0, 1];R''), the functionals and K"' are respectively regarded as M'' x M''-valued and 
M'^ (g) M'' (g) ^''-valued L^-functions on [0, 1]: 



if" = (a-(s)l; 



0<s<l. i=l,...,d, a=l r 



and 



(^"(s)l^; 



0<,s<l, i=l, 



, a,p=i,...y 



The sequence {ij}j=o,i,... j" £ N) is a triangular array of numbers such that tj — depending on n and 
that = to < ^1 < • • • < tj" = 1- Functional /f"(s, r) is defined as the R'' (g) R' (g) R'-valucd function 



ii'"(s,r) 



for (s, r) e [0, 1]^. We write 



i = l,.-.,d 
<,fJ=l,.-.,r 



Suppose that K" and /C" are progressively measurable. More strongly, we assume the strong predictability 
condition that K'^{s) is J^t^_j -measurable for s G {tj-i,tj]. Corresponding to the kernel if", we consider M" 
given by 



K-{syJJ2 k-{ry^pdwf!Adw^ 



(11) 



i=l,...,d 



Write C" = {M^)t and C„ = (M")i. Suppose maxj \Ij\ — o(r„), where Ij ~ {tj^i,tj], and that the sequence 



of measures 



(12) 



weakly for some measure fj, on [0, 1] with a bounded derivative. We will assume that J2j l-^jl^ — 0{1). In 
this case, 



•^"(^)a^(E /' K-{ry^^dw^^) ds] 



l^dw!^ 



and the in-p limit of C„ will be 



Co 



i5]/'i?-(t,t)-,if-(t,i)L^,MdO 



ii ,22 — 1,. ..cZ 



K°°it,t)(g) K°°{t,t) ii{dt) 



(13) 



under the conditions we will assume, where K°° is the limit of if", and Tr* is the trace on (R''(gR'')*g)(R''g)R''). 



Let q G (1/3,1/2). Let A = {(s,r); < r < s < 1} and A" = Uj{(s,r); tj_i < s < t < tj}. 



Let 
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[Al] (i) K" £ 3i+i{H (E) H (E) M.'^) and a representation density of each derivative admits 



ess. sup 

ri,...,rfc,e(0,l), 
(s,r)eA", TiEN 



Dr 



R^{s,r) 



< CO 



for every p € (1, cxd) and k < £ + 1. 
(ii) For every 77 > and p G (1, 00), 



sup 

se(o,i),eeS'' 







r i4e®2] 1 


-1 


1 




(1 - s)l+'7_ 





< 00. 



(iii) * For every p > 1, 



and 



sup sup |lif"(s,r)-i^"(t,_i,t,_i)|l,,p = Oirl") 



j (a,r)eA" 



supsup||i?"(i,_i,t,_i)-i?°°(t,t)||,,p = 0(r29) 
j t&ij 

as n — >■ cx). 

Remark 5. Condition [Al](iii) is [Al](iv) of [8]. In [8], [Al]'' was "Conditions in [Al] except (iii) hold." 
Remark 6. (i) In typical cases r„ = n^^/^ so that r^^'^ = nr^ > n^i for n > 1. 
(ii) Under [Al] (i), for every p e (1, 00) and k < £, 



ess. sup 

ri,...,r^,se{OS) 



\Dr„...,rJs\ 



1~S 



< 00. 



(iii) As for [Al] (ii), we need the nondegeneracy of the derivative of in s, or a large deviation argument, in 
order to control eTip{^{C^ - C^)[u'^'^])) for s near 1. 

Let Coo e -F(ri;R'^ ® R'*), G ^(17; R'^) and Foo G J^{n;R^^). Set C„= r-i(C„ - Coo), Wn^ r-^Wn - 

o 

Woci) and Fn— {Fn — Foo). Consider an extention 



000 



l'^ such that 



of in,T,P) by a probability space (rj, J",F). Suppose that G J"(rj; C([0, 1]; R'')), A^oo G J"(r2;R'^), CooG 

J"(fi;M''«)M''), ^ooG -7^(f^;K'') and Fo^e Tin^W^^). 
For JF = F y (j[M^], there exists a measurable 

Coo(w,Afoo) = £^[Coo|-^]. 

Similarly we define Woo(w, z), Foo{ijJ, z) and Nao{^, z) by 



W^oo(a;,Moo) 

i^oo(a;,Moo) 
7Voo(a;,Moo) 



E[Woo 1-^] 

^^[^'oo 1-^] 
^^[A^ool-^]. 



Further, we introduce the notation 



Coo (z) = Coo (1^,2) := Coo(w,2; - VKoo) 

and similarly Wooi'^^z), i^oo('^,-z) and Noo{uj,z). 

Let Sn be a positive functional defined on O. It is said that a functional c : D, X R'^ ^ M is a 00 -polynomial 
if c is a polynomial of 2; G R'' with coefficients in ©s.oo- 



''Condition [Al] (iv) of [8]. 
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[A2] (i) Foe € Df+i.^(M''i) and Woo e Be+i,oo{M.'^) ■ 

(ii) Fn e D£+i,oo(M'^0, Wn G Df+i,^(R''), Nn e Df+i,oo(M'^) and s„ e Df,oo(M). Moreover, 

sup ||| Wn \\i+l,p + II F„ \\£+l,p + ||A'"„||£+l,p + ||s„||f,p| < oo. 

for every p > 2. 

(iii) C^, W^, {j,k = l,...,d) and Pj;^ {I = l,...,di) are D^o^oo-polynomials in z G R"^, where 
4 = 2[(rfi+3)/2]. 

(iv) (A/.",iV„,(7„,Ty„,F„) (M.°°,7V^,Coo,T^oo,Foo). 

(v) ForG-T^oo andFo,, 

ess. sup ||^ri,...,rfcG||p < OO 

ri,...,rfee(0,l) 

for every p e [2, oo) and k < £ + 1. Moreover, r i— >■ D^G and (r, s) i— > Dr^sG are continuous a.e. with 
respect to the Lebesgue measures. 



The nondegeneracy of (M" + Woo,Frx,) will be necessary. 
[A3] (i) There exist a sequence (t„)„gN in [0, 1] with sup„ i„ < 1 such that 

suptx^ P[det cr(^/n_i_vi/_^ j?^) < Sn] — 0{r'^) as n oo for some v > £/3. 

(ii) For every p> 2, limsup„^oo E[s:^p] < oo. 

Remark 7. The nondegeneracy det G^^^ G C^p>2L^ follows from [Al] (ii). Indeed, it implies 

sup P[Goo[e^'] <e] < sup WMe'^XX ^ ^p^" > «) 

for some constant Gp for every p > {]. Then the desired inequality is obtained; see e.g. Lemma 2.3.1 of Nualart 
[4]. 

Now we recall the martingale expansion ([7, 8]). For Z„ given by (10) and (11), the random symbols are 
specified as follows. The adaptive random symbol is 



a{z,iu,iv) = ^Coo{z)[{iu)^'^] + Woo{z)[iu]+ Noo{z)[iu]+ Fooiz)[iv]. 
For the double stochastic integral in question, the anticipative random symbol is given by 



(14) 



cr(m, iv) 



1. 



-Tr* 



K°° {t, t) [iu] ® at,t [iu, iv) fj.{dt) , 
where the random symbol crt,f,(iM, iv) has the expression 

atAiu,iv) = (iDtWooM - ^AGoo[m^'] + zAPooH) ® (zAW^ooM - ^AGooK'] + iA^ooH) 

+ {iDtDtWoo[u] - ^DtDtCoo[u^^] + zAAPooH) 

with the representation densities of the Malliavin derivatives of functionals. The approximate density pn{z,x) 
is defined by 



(2;VKoo,G„o)4(-Foo) 



a{z,dz,dx)*\ (l)iz;Woo,Coo)SxiFoo) 



Pniz,x) = E 

for the full random symbol 

a ~ a + a. 

Note that u is d-dimensional here. Class £{M,j) will be abused for functions on M"*. 



(15) 
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Theorem 2. Suppose that Conditions [Al] , [A2] and [A3] are fulfilled. Then for any positive numbers M and 
1, 



sup 

/e£(M,7) 



E[f(Zn.,Fn)\ - / f{z,x)pn{z,x)dzdx 



o(r„) 



as n —)• oo. 



4 Some details of derivation of the expansion for the quadratic form 

In this section, we will give somewhat detailed proof of Theorem 1, which was originally presented in [8]. 
4.1 Stochastic expansion 

We will work with the Ito stochastic integral equation (1). The following lemma gives a stochastic expansion of 
the targeted variable Z„ of (3). 

Lemma 1. Z„ has the following stochastic expansion: 
where 



j—l Jtj-l/\t Jtj-1 



and 



Nn = Qn^ct-_^^at^_^a\f^_^ / / / dwudwsdwt 
j—i '^tj-i Jtj-i Jtj-i 



+2'^ct^_,bt^_,at^_, dwt + 2n^cf^_icrt^_,cr|^^l / {t-tj.i)dwt 
j=i Jti-i j—i Jt. 

n n 

j — l J tj-l J tj-l 

1 [0] 2 1 [1] [1] , /i\ 



2n^ ^.-1 t.-i n 

Here Om(1) denotes a term of o{\) as n ^ oo with respect to Jigp-norms of any order. The families {Af"}tg[o ij 
and {Nn}neN o,f^ bounded in every IJ>s^p-norm. 

It is possible to obtain the above lemma by somewhat long computations. See Section 5. 

4.2 Random symbols 
4.2.1 Adaptive random symbol 

The discrete filtration F" ~ (-^")te[o.i] with .F" = J^[nt]/n will be used. The predictable quadratic covariation 
for F" is denoted by ((•, •)). Note that ((•, •)) depends on n. Let Hi{x) = x and H2{x) = {x^ — l)/^/2. Denote 
A.jW — wtj — wtj_-^ , which depends on n as well as j. The discrete version of M" is given by 



Af,2'" = ^ E V2a{Xt^_,)H2{V^A,w). 
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For C'i ■■= V^(Ct" - Cf=) , we have 



) — (s — ij-i) (ds 



) — (s — ij-i) /ds + Op I 



Here the supremum of "Op(?i-i/2)" in < e [0,1] is of Op(n-i/2) ^jth respect to LP -norms. Therefore, the 

o 

principal part of C" is F"-inartingale 

" = 51 / ' 4ny^a(Xt^-_j2j f dwr)" - {s - \ds 



The discrete version of w is denoted by — W[„t]/„. 
By a similar argument, we have 

~ ~ o 

In the present situation, Woo{z) = and -^00(2) = 0. We need to identify the limit (Moo, Cooj-^00) to write 
the adaptive random symbol. The "martingale part" of Nn with respect to F" is given by 

r*3 /■* 

TV" = 6n ^ ct,_, crf,_, ci^' , III dwudwsdw. 



j:tj<t Jt]-! Jtj-i Jtj^i 

/tj /"Ij 
dwt + 2n Ct^_iO-t,_iCr|^^Li / (t-tj^i)dwt 

■ti<t •'tj-i •'tj-i 



-n 

j:tj<t 
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Then 



((M«'",M«'")) 



((u;",iV))t I hds 



T2,n Arnv 





M 

n 




2 




71 




16 




in 

3 




/ ' 

Jo 






o 
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3ri 




0, 


-^p 


0, 


->p 


/' 



((iV",iV"))i / gfrfs 

Jo 

as n — > cxD for each t € [0, 1], where R+-vahied process qt takes the form 

% = P(ct,c[',bt,at,crl') 

for some polynomial p; it is possible to give an explicit expression of p, however we do not need the precise form 
of qt later. The orthogonality of Af^'", M^'" and to any bounded martingale orthogonal to w is obvious, 

o 

thus with the representation of C", Af" and Nn, we obtain 

(Afoc,Coo,iVoo) (/ V2a{Xs)dBsJ ^a{XsfdB,+ ( \a{X,fdB'^, 

^Jo Jo -J Jo "J 



/ ksdWs+ f ^M^dB'^ + f hsdi 

Jo Jo Jo 



where {B, B' , B") is a three-dimensional standard Wiener process, independent of J^, defined on the extension 
n. Since 



Noo{z) = f ktd'Wt+ / htdt, 
Jo Jo 

the random symbol a{z,iu,iv) is given by (5). Moreover we see Condition [A2] holds. 
4.2.2 Anticipative random symbol 

Let us find the anticipative random symbol a(iu, iv). Recall that a{x) — a{x)'^, 

Cr = 2 / a{Xt)dt, Coo = 2 / a{Xt)dt, ^ f f3{Xt)dt and W^, = 0. 

Jo Jo Jo 

To describe <Ts, s{iu,iv) more precisely, we consider the random symbol fjs.r{iu,iv) that admits the expression 

(Ts^r{iu, iv) 

= u^^ a'{Xt)DrXtdt(^-u' a'{Xt)DsXtdt + i {Xt)[v]D sXtdt) 

+ (-^2^ a'{Xt)DrXtdt + t P'{Xt)[v]DrXtdt)(^~u^ a' {Xt)DsXtdt + i j ^ 13' {Xt)[v]DsXtdt) 
+ f {a"{Xt)DrXtD,Xt+a'{Xt)DrD,Xt}dt + i f {(3" {Xt)[v]DrXtD,Xt + (3' {Xt)[v]DrD,Xt}dt^ 
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for r < s, where the prime ' stands for the derivative in xi € M. The random symbol as^s is the hmit 
Irnir^s o's.r{iu,iv), that is, (7), and this gives the anticipative random symbol (6). 

4.3 Nondegeneracy 

Here we will consider the nondegeneracy of (M", Foo) in Malliavin's sense. Let 

r]j{t) = ^/n{w{tj At) - w{tj_i A <)) 

and 



Then 



Ml' 



^ 3 = 1 

The representing density of the Malliavin derivative of is 

n 

1 " 

=: Di{n,t)r+D2{n,t)r. 



<t} 



Now 



1 " 

D2{n,t)r - — ^a'(Xt^.„jA.Xt^_iOWl{r<t,-i<t} 

n — 1 

= 0-'iXte)Ce+lit)DrXtillr<te<t} 
£=1 

n-1 £ 

V fcl j = l 
-|^n— l^n — 1 \ 

= ^ E ( E ^'iXt,)^,+,{t)DrXt, ) l,^.(,)(r) 



where Ij{t) — (ij-i A t,tj A t]. Hence the Malliavin covariance of M" is 

ftjM - 



for li G M, where we read ^ 



fc=n+l 



0. Since 



p'{Xs)DrXsds, 
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we have 



ai2{n,t)[v] {M^,F^[v])h 



nlj /\Z n I 1 



P'{X,)DrXsdsdr[v\ 



for w e . Let 



and let 



®2^ _ 



(j{n, t) 



Let 



0-11 



1 " 



0-11 ai2{n,t)* 
(712 cr22(i) 



- 2 n /.(jAt 



dr. 



(ir, 



'7i2(«,i)H = 



J r=t,^, At J s=r y/n 



and 



(5-ii(n,t) a-i2(n,t)* 
ai2{n,t) (722 (t) 



We shall show 



|(7(n,i) = 0(n 2) 



(16) 



for every p > 1 and t G [0, 1] (in particular, for t = 1/2). For this purpose, we need a lemma. Let X denote the 
set of sequences j'"^-* = [JIi ]) of multiple Ito stochastic integrals taking the form 



si 



where {a^n\ * = j = ^,---,n, n G N} is a family of progressively measurable processes, bounded in 

]D)°° = r\s,p^s,p- In the following lemma, Jn,]l ■ ■ ■ Jnjl^ j'^kl ' ' ' j'nkl ^"^^ ^^'^ ea,ch. of them has Onj^i 

which may possibly differ from those of other indices i^'s and /i's even if the values of indices coincide each other. 

Lemma 2. Suppose that 



sup 

je{i,...,Ti}, nEN, 7e{o.i,...,m}. 

T-i,....,r^,se[0,l] 



for all O = a|*^, l}"*^- and a^*^- ^{s), and for every p > 1.^ Then 

(a) Suppose that b^^l, are Tt^_j^ -measurable. Then for 1^1, ...,i/m, fJ-i, fJ-q & 



1 

— fo^''' /'^'■' 

/ ^ n,kq n.kq 



^7 = denotes the case with no derivative. 
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(b) For iyi,...,iy„, € 



1 " 

± y E 



= o 



1/2 



T/ie constants in the above estimates depend only on the given supremums. 
By using smoothness of appearing functions, apply Lemma 2 (a) to 

rtjAt 

bi% = n a\Xt,_,) / DrXt,_,dr {j < k) and J^l = ^-(0 



to obtain 



te[o,i] 

Furthermore, Lemma 2 (b) apphed to 



sup crii(r;-, i) — (Tii(n, 



element of 2a{Xt._-^) x n 



tj At 



r—tj_i/\t J s—r 



p'{Xs)DrX,dsdr[v] and J^^f = 7^j{t) 



yields 



sup ||o-i2(n,t) — (Ti2(n,i)| 
te[o,i] 



as n — > 00 for every p > 1. Consequently, we obtain (16) 



Let m„ = 71 ^ '7j(l/2)^- For a positive number ci, define s„ by 



- det 
2 



where ^ K ^ [0, 1] is a smooth function such that ip{x) = 1 if |a::| < 1/2 and "ipix) = if |a;| > 1. Then 
Sn > if m„ < ci, and s„ > 2^^ det(7(n, 1/2) otherwise. Thus, it suffices to show 



sup£;[l{„„>ci}(detCT(n, 1/2)) ^] < 00 (p > 1) 



(17) 



for the nondegeneracy 



sup i?[s„^] < 00 



(18) 



for every p > 1. Following precisely, e.g., the proof of Lemma 2.3.1 of [4], in order to obtain (17), it is sufficient 
to show that for every p > 1, there exists a constant Cp such that 



sup P[mn > ci, a{n, 1/2) [u®^] < g] < CpC^ 

UGK1+<*1 :|u| = l 



(19) 



for all e £ (0,1) and all n G N. [The reasoning there is valid even for the measures l{-„j^>cj].(iP in place of 
P.] Here we use L°°^-boundedness of {a(ri, l/2)}„gN. Besides, for a while we shall assume the nondegeneracy 



condition: for some constant Cr, 



sup P[a22(l/2)[«®2] <e] < CpcP 



(20) 
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for all e €E (0, 1). Suppose that |u| = 1 for u = {u,v) € 
0-22 for (722(1/2). Let p > 1. When \v\ < e^/^ 



For simplicity, we write (jy for aij{n, 1/2) and 



P[m„>ci, a(7i,l/2)[u®2] <e] 

P[m„ > ci, CTiiu^ + 2(5-12 M + cr22b®^] < e] 



< P 



l„ > Ci, 2|(712U| > 



2ei/8 



P[to„ > ci, (Tiiu^ < 2e] 



< ^'[4|ai2| > e~^/^^|M|] + P[m„ > ci, cthu^ < e^/^^] + P[m„ > ci, anu^ < 2e] 



< 



(21) 



uniformly in ri G N and u satisfying |u| = 1 and |ti| < e^/^, where C2 '■— inlx > 0, since (1 — e^^^y^^ < 

\u\ < 1 and 



for any eg (0, 2 ^) on the event {m„ > Ci}. 
We will assume \v\ > e^/^. From 



aiiu^ + 25-12H + cr22[f®^] > <5-ii^|a-iiCT22[w®^] - (cti2[u])^|, 



it follows that 

P[m„>ci, a(n,l/2)[u®2] <e] 
P[to„ > Ci, CTuU^ + 2CT12M + (722 b®^] < e] 
< P[au > +PK > ci, ^ii(T22b®'] - ('^12H)' < e'/']. 

On the other hand, we have 



®2i 



('5-12 H)' 



1 " 

-E 



2 



j=l •'tj-iM 

E 



^ a'(X,,_ja(i)^rXt, 



dr 



p'{X,)DrXsds[v 



> 



-E 



2aiXt^_,)r,,it) 



k=j+i 
2 



dr 



P\Xs)DrXsds[v] 



dr 



for i = 1/2, where we used the Schwarz inequality. Hence, 

ana22b®']- ((^i2H)' > Ac.cl cT22il/2)[v'»^] 
on {m„ > ci}. Combining (20) with scaling v t-^ e~^/^w and (23), we obtain 

sup PK > ci, a^ia22[v^'] ~ (^i2H)' < e'/^ 

„gRdl:ei/8<|i,|<i 

< sup P[4cic2a22(l/2)[t;®2] <e3/4] 

t,gR<ll:£l/S<|i,|<i 

< sup P[a22(l/2)[«®2] < i4ciclr\^^^] 

veR''l-:\v\ = l 



(22) 



(23) 



(24) 
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for every p > 1. By connecting (24) to (22), we obtain 



sup P[m„ > ci, o-(n,l/2)[u®2] < e] < e 

uGRi+<il: |u| = l, ei/8<|t,|<l 



P 



(25) 



Thus from (21) and (25), we obtain (19) and hence (18) under the assumption (20). 

We consider the (1 + di )-dimensional process X = (Xt)tg[g satisfying the stochastic differential equation 



dXt 



Vo{Xt)dt + Vi{Xt) o dwt, Xo = (Xo, 0) 



in the Stratonovich form. Then the Hormander condition [H2] together with the compactness of suppP''''" 
ensures that for t € (0, 1] and for every p > 1, there exists a constant Cp such that 



sup P 

vGK1+'*1 :|v| = l 



D,Xi DsXidsW'^^] < € 



< CpeP 



(26) 



for all e G (0,1). See Kusuoka and Stroock [2, 3], Ikeda and Watanabe [1], Nualart [4]. In particular, (20) 
follows from (26) applied to v = (0, v) for v G K''^ . 

We choose ci such that 2ci < 1/2 = limn^oo E[mn]. Since 



by definition, we see that for every K > 0, 



(Ti2{n,t) CT22(1) 



> cr(n,t) 



sup P[detai^Mi',F^) < Sn] < P[det cr(n, 1/2) < s„] + 0{n ^) 

t>l/2 

< P[deta(n,l/2) < 2s„] +0(n--^) 

< P[m„ <2ci] +0{n-'^) 



(27) 



as n — > oo. Here we used (16). Properties (27) and (18) verify [A3]. 
4.4 Proof of Theorem 1 

It is easy to verify Condition [Al] under [HI]. Conditions [A2] and [A3] have been proved in the preceding 
sections. Now we can apply Theorem 2 to obtain Theorem 1. 



5 Proof of Lemma 1 



Let 



We have 



/(i.i) = = a' a, f^'^^ = a™ = a'b + ^a"a\ 



3 



(28) 
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In the decomposition (28), 

- ^E%-.(/ej^ + E2^*.-./e^./£?r f d^sdt 



j j 



as n — ^ oo. Next, 



= ^Y.'^t.^-fhUii r f dwsdwt+2Y,ct,.jt'^fiii r f r dw^dw^dwt 
= 2Ect,-,/sji;^, r r dw,dwt+6Y.''t,.jt'^fiii r f r dw^dw^wt 

j J tj-i J tj-i j J tj-i J tj-i J tj^i 

+2j2^t,^Jt-ljtlL r f dsdwt + 2j2c,^_J^l:';>f^ll P f dsdwt+oJ-\. (29) 



16 



J tj-.i J tj-i j J tj-i J tj-i J tj-i 



= (/Si)'+r 2/«/^^w.+/ 2/w/^^)ds+/ 



t,_l 
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Now 

/ Cta^dt — / cta^dt 
j Jtj-l Jtj-l 

= ^E^*.-.(iSJ^+E2c.-.iS^.ie:V' f "+;!^E%-./ff.ifc^^ + iE^^.-.(*:^)^ 
j Jtj-i Jtj-i 

4 E + i e (^e + f ^) • (31) 
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From (28)-(31) we obtain 

J Jo 



Vudwgdwt 



„2 ^'^ti-iftjljtjli + E"^*J-i('^*j-i^) 

+2y.c,^_j^ij:i p f d^.d^.+6^c,^_ji::^f^iir r r d^. 

j J tj-i J tj-i j J tj-i J tj^i J tj-i 

+2E^*.-if r r dsd:w,+2Y^c,^_ji:^fiii p f dsdw 

3 3 

- { ^ E (if + E 2c.. i^. it? /' f '^-^t + -i E c... /t . + i E ) 

3 3 i 3 

+E^-.(/tj'r r '^-^^^ 

— Vr(l) f(l) f (I'l) , J_Vr(2) rf(l) ^2^1 -^1 



f(2) n2 



— ' 'J^' 2^2^ — ' ■■'-^ -J-^ J ■ \n 

i 3 

L J J tj-i J tj-i J tj-i 

+^T.<^^.-J^-JSI r f dsd^, + 2^2 c,^_J^:^ fill p f dsdn., 

-E<'-.(/tj'r r '^-^'^^ 

2^2 Z^%-lWtj-i/' ^2 ^-iJtj-i-'tj-i J 



o„ I - 
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Thus we obtained 

Ut 



V ^ Jo 

J t-i -^T J t-i ^-1 



dwudwsdwt 

tj — 1 " ij — 1 ij ~ 1 



-\-2n^Ct._^bt._^atj_^ / / dsdwt 2n^Ct^_^at^_^(7^^^_^ / / dsdwt 
j J tj-l J tj-l j J tj-l J tj-l 

1 1 

E^-i4-i t f '^^^d^ 



-n 

j 



2 

J 3 

dwgdwt 



ftj rt ru 

dwudwgdwt 



j tj-l J tj-l 

+'^^^t,^ibt^_,at^_, I dwt+2n'^ct^_^at^_^^(^t}^^ I I dsdwt 

j j J tj-l J tj-l 

'3 j 

-'^E^-i'^t.-i / / '^^^^^ 

1 [0] 2 1 [1] [1] , n\\ 

Moreover, we can obtain Bs^p-estimates of the residual term as weU as each term on the right-hand side. 
This completes the proof of Lemma 1. 
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